Controlled Nanoparticle Formation by Diffusion Limited Coalescence 
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Polymeric nanoparticles (NPs) have a great application potential in science and technology. Their functional- 
ity strongly depends on their size. We present a theory for the size of NPs formed by precipitation of polymers 
into a bad solvent in the presence of a stabilizing surfactant. The analytical theory is based upon diffusion- 
limited coalescence kinetics of the polymers. 

Two relevant time scales, a mixing and a coalescence time, are identified and their ratio is shown to determine 
the final NP diameter. The size is found to scale in a universal manner and is predominantly sensitive to 
the mixing time and the polymer concentration if the surfactant concentration is sufficiently high. The model 
predictions are in good agreement with experimental data. Hence the theory provides a solid framework for 
tailoring nanoparticles with a priori determined size. 



Polymeric nanoparticles (NP) are gaining increasing atten- 
tion because of their numerous applications in, for instance, 
physics, chemistry and medicine lM- The NP size and size 
distribution are the key parameters often determining their 
functionality. Therefore, one of the main experimental chal- 
lenges is to prepare NPs with well controlled dimensions 
tuned for a particular application. Models of NP formation, 
allowing one to steer the NP preparation process in the right 
direction, would simplify the size control significantly. 

A high level of control over particle size is required in, for 
example, targeted delivery (e.g., oncology). Size influences 
the circulating half life time and is crucial for selective cel- 
lular uptake: NPs between 50 and 200 nm in size are de- 
sired in passive cancer tumor targeting as they are too large 
to harm healthy cells but small enough to penetrate into the 
diseased ones. In brain imaging, fluorescent dye loaded par- 
ticles of about 100 nm with biocompatible polymer coatings 
are used because they produce small, sharply defined injection 
sites and show no toxicity in vivo or in vitro flj-Q]. 

Although there are several methods for NP preparation, 
only few of them permit high level of control on the parti- 
cle size and the particle size distribution J2]. Often, a water 
insoluble moiety (e.g., a drug or a dye), needs to be encapsu- 
lated into a carrier polymer and protected by an emulsifying 
agent, which also makes the NP water soluble. In particular, 
the so-called nanoprecipitation method permits preparation of 
nearly monodisperse NPs in a very simple and reproducible 
way |5j|. Typically, an organic phase, which is usually a dilute 
polymer solution, e.g., PCL in acetone, plus the hydrophobic 
moiety to be encapsulated, e.g., a drug or a fluorescent dye, 
is injected by pressure into an aqueous solution of the emul- 
sifying agent, Figure [Tf a )- As the organic solvent is chosen 
to be water-miscible, rapid quenching (towards poor solvent 
conditions) of the hydrophobic polymer and the drug in water 
takes place. This results in coalescence of the polymer and the 
drug into submicron particles decorated by surfactant J2,0]. 
Alternatively, a block copolymer can be used to combine the 
polymeric drug carrier and surfactant roles into a single com- 



ponent. To promote a faster and better controlled mixing, an 
impinging jets mixer, Figure [Tib), has been employed in this 
case Jj,|7D. 

Despite the simplicity of the experimental method, one still 
lacks a comprehensive theoretical model that allows to pre- 
dict how particle size depends on the materials and process 
parameters, in particular on concentration. Therefore in many 
practical situations investigators still resort to simple empir- 
ical correlations [8|] or to statistical methods such as experi- 
mental design 06[]. More advanced theoretical methods, such 
as Brownian dynamics simulations [9] or population balance 
coupled to CFD simulation ifioll . do provide very valuable in- 
sights into the kinetics of mixing and rapid assembly upon 
quenching but do not permit formulation of a simple yet phys- 
ically meaningful analytical relationship between the experi- 
mentally relevant parameters and the NP size. Such a rela- 
tionship would be extremely useful in designing NPs with a 
priori determined size as it would allow one to avoid a very 
laborious trial and error process. 

In this Letter we formulate a model of the nanoprecipitation 
process. The model grasps all essential features of the process 
and, at the same time, provides a simple analytical expression 
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FIG. 1. Scheme of a pressure driven injection device used in @| (a) 
and an impinging jets mixer used in 17] (b). Fluid A is the organic 
phase comprising solvent, carrier polymer and the drug, fluid B is an 
aqueous solution of emulsifying agent. 
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for the NP size as a function of the mixing intensity and the 
surfactant and polymer properties. 

We restrict ourselves to a bi-component system: a dilute 
solution of a hydrophobic polymer is injected into a wa- 
ter/surfactant solution. Because the solvent and water are cho- 
sen to be well miscible, a rapid quench of the hydrophobic 
polymer in water takes place and the polymer particles start to 
coalesce upon encounter to form larger particles. In parallel, 
the surfactant molecules, also subject to Brownian motion, ad- 
sorb on the surface on the newly formed polymeric particles 
and make their coalescence progressively more difficult until a 
stable situation is reached. As the NPs formed represent a sys- 
tem in a kinetically frozen state, their parameters will depend 
strongly on the system kinetics, which includes at least three 
processes, namely (i) mixing of the polymer plus solvent with 
the aqueous surfactant solution taking place on the time scale 
Tmix', (ii) coalescence of the hydrophobic polymer particles in 
a hostile water environment, characterized by a time t c i s ; and 
(iii) protection of the polymeric NPs by the surfactant taking 
place on the time scale r pro and bringing the system into a 
kinetically frozen state. 

Let us first consider the limit of 'very fast mixing' r m j x — » 0, 
in the absence of surfactant. The system then initially consists 
of collapsed polymer molecules homogeneously distributed in 
water. These molecules will diffuse, collide, and stick. If they 
would be hard particles, this would lead to fractal aggregates, 
for which well-known growth laws have been developed. This 
case is commonly known as 'diffusion limited aggregation' 
(DLA). However, as our particles are liquid like, they will co- 
alesce to homogeneous spherical particles rather than forming 
fractal aggregates, so that we deal with 'diffusion limited coa- 
lescence' (DLC) lllll . For this we have Smoluchowski theory 
1 1211 with a rate K = AnD'R' ', where D' and R' are the sum of 
the diffusion coefficients and the radii of the reacting species, 
respectively. Hence, in a mean field approximation, the poly- 
mer particle concentration c p is governed by a simple equation 
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where the Stokes-Einstein expression D p = k B T / \6nrjR p ) for 
the diffusion coefficient of a particle in a fluid with viscosity 
i] has been used. The factor h accounts for the probability that 
a collision leads to a coalescence event. 

The surfactant adsorbed on the particle surface influences h, 
as it reduces the probability of a coalescence event to occur. 
Hence, h is a function of the fraction of the particle surface 
protected by the surfactant, h = h\n(t)a 2 /(4nR 2 (t)y}, where 
«(f) denotes the average number of surfactant molecules ad- 
sorbed on a polymer particle with radius R p at time f, each 
surfactant molecule covering a surface area a 2 . 

Since we are dealing with coalescence rather than aggrega- 
tion, there is a direct relation between particle mass and parti- 
cle radius R p leading to the mass conservation law in the form 



cM)Rt(t) 



c p oi?pQ. Here, Rpo and c„o are the size and the 
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number concentration of the polymer particles immediately 
after the mixing took place. 



In the absence of surfactant h(-) = 1 and ([U can be recast 
in terms of the particle size yielding Rp(t) = R pQ (1 + f/r c i s ), 
with an encounter and coalescence time 



(2) 



Assuming strong favorable interaction between polymer 
and surfactant, the surfactant-polymer coagulation can be 
treated in a similar manner 

dc* ee 2k B T I 1 1 \ , x r 

where R s and cjj ree are the diffusion radius and the concen- 
tration of the free (not adsorbed) surfactant molecules; h s - 
h s (na 2 /(4jtR 2 )) denotes the probability of adsorption. Equa- 
tion (f3j is a straightforward generalization of ([TJ, where a rel- 
ative diffusivity has been introduced as a sum of the polymer 
particles and the surfactant molecules diffusivities (see lfl3ll 
for more details). Also, the reaction radius is assumed to equal 
R p + R s . 

An interesting observation at this point is that T pl0 ~ t c i s 
and, hence, collision rate of the polymer particles and their 
protection by the surfactant go at approximately the same 
pace. Note, that we have neglected surfactant micellization 
by assuming that the surfactant molecules bound in surfac- 
tant micelles behave similarly to the dissolved ones, at least in 
what concerns their agglomeration with polymeric NPs. 

To describe the kinetics of coagulation, the exact functional 
form of h and h s must be specified and the equations ([TJ and 
(O need to be solved together. In fact, as we are only inter- 
ested in the final particle size and not its time dependence, 
we can divide ([TJ by @ yielding a single differential equa- 
tion for c p as a function of cjj ree . Computing the exact form of 
h(-) can be quite involved, although it is clear that h(Q) 1 and 
h(l) si 0. The same holds for h s . To simplify the matter signif- 
icantly, we assume the surfactant adsorption not to influence 
the coalescence of particles until the particles are saturated 
with the surfactant and the coalescence is stopped completely 
lfl4fl and take h(0 < x < 1) = 1 and h(x > 1) = and the same 
for h s . Such a choice does not change the scaling of all the 
important quantities but implies that coagulation process stops 
when n = AnR 2 p la z and, hence, c f ^ nd = c s0 -4n{R p nA ) 2 cp nA la 2 . 
Solving the differential equation for c p as a function of cf ee 
explicitly and making use of the above relation between the 
end values of the concentrations and the radius, one derives a 
transcendent equation for the ratio £ = R p nd /R p o between the 
final and the initial particle size 
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where a — R p q/R s is the ratio between the initial polymer 
particle size and the diffusion radius of the surfactant and 
k = 47rR 2 pQ c p o/(a 2 c s o) is the ratio of the total initial surface 
area on the polymer particles to the maximum area surfactant 
molecules can occupy and block. 
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Analytical solutions of are found for the limiting cases 
of an excess of surfactant, k <k 1, and if surfactant is scarce, 
k » 1: 

I K 

fl™ d = /? p0 x ^ 3/4 + a + a- 1 (5) 
k if k » 1 . 

This leads to a simple interpolation /?™ d ^ /? p o(l + *)> which 
is surprisingly close to the exact numerical result. 

A very peculiar implication of the fact that r c i s ~ T pro , as 
pointed above, is that the final NP size does not depend on 
the mobility of polymer or surfactant molecules. The only 
dominating factor in the 'fast mixing' limit, when T m ; x <K r c i s , 
is the surfactant concentration. 

Let us now consider the other limit, r m i x > t c i s , which is ap- 
parently characterized by a very fast particle aggregation on 
the time scale shorter than the typical mixing time followed 
by stabilization of the NPs' size at the times t > T m ; x . In- 
deed, at the very beginning of the process, the polymers are 
present as isolated chains in a good solvent. As the solvent 
quality drops the polymers instantaneously collapse. Subse- 
quent collision of collapsed chains leads to coalescence fol- 
lowing the kinetics prescribed by ((TJ with h = 1. Hence, 
the particle size at the end of the mixing, t ^ r m ; x , reads 
-Kmix - Rpo (1 + Tmix/Tds) 1 ^ 3 . At longer times, there is enough 
time for the surfactant to adsorb onto the surface of the coa- 
lescing polymer-rich dropets. Then the system finds itself in 
a well mixed state and its kinetics obeys the set of equations 
([TJ and (f3]) as discussed above, but R m [ x must be used as the 
'initial' particle size in ©. This two-step process leads to a 
final expression for the polymer particle radius in a kinetically 
frozen state 

/ T ■ \' /3 

^=R p0 a+K)\i + ^\ . (6) 

It is characterized by a plateau at small r m i x /T c i s , where the 
NP diameter is independent of mixing or encounter and co- 
alescence time and is totally governed by the surfactant con- 
centration (parameter k) with the smallest particles obtained in 
excess of surfactant. The other regime, T m i x /r c i s > 1, shows a 
typical 1/3 power law behavior and is dominated by the mix- 
ing efficiency. 

Based on © it follows that for typical experimental con- 
ditions, i.e. excess of surfactant and relatively slow mixing, 
the final NP size depends mainly on the mixing time and the 
initial polymer mass concentration c mp , R p nd oc (cmpTmix) 1 ' 3 , 
and is independent of the polymer molar mass. Only a mi- 
nor dependence on the molar mass of emulsifying agent can 
be observed indirectly via T m ; x , which can be sensitive to the 
viscosity of the surfactant solution. The same holds for the 
temperature. 

To appreciate the formula §6^ we compare its scaling pre- 
dictions to our own experiments as well as to the data avail- 
able in the literature. In the nanoprecipitation experiments 
performed in our lab, PCL (CAPA 6250 supplied by Solvay) 
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FIG. 2. Size of PCL (molar mass M„, = 25 kg/mol) NPs prepared 
from an acetone solution quenched in an aqueous Pluronic solution 
as a function of the initial polymer concentration. Solid line shows 
a fit to the data. Also the data from Molpeceres et al. [ 6] are shown. 
Inset: Size of NP vs molar mass of PCL, at concentration c mP ct = 
5 mg/ml. 
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FIG. 3. Master curve of the size of the core of the diblock copolymer 
NPs versus the rescaled mixing time, following the scaling predicted 
by l|6}. Inset: The original data from Johnson and Prud'homme [21]. 



has been used as a carrier polymer and Pluronic (PF127 sup- 
plied by BASF) as a surfactant. PCL solution in acetone were 
quenched in a 1 wt% PF127 aqueous solution with a device 
similar to the one depicted in Figure [TJa). The hydrodynamic 
particle diameter D/, has been measured by dynamic light scat- 
tering. 

As can be seen from Figure |2l our results compare favor- 
ably to the data available in the literature |0] for the same 
system. As the experiments are performed in the T m ; x > r c i s 
regime, the scaling obeys the 1/3 power law as expected. To 
check the molar mass sensitivity, additional experiments have 
been performed where PCL molar mass has been varied be- 
tween 2 and 80 kg/mol. The diameter was, however, hardly 
affected by the molar mass [see inset in Figure [2), in accord 
with the theoretical predictions. 

The data presented only cover the T m ; x > r c i s regime and 
neither reach a particle size saturation limit at the very fast 
mixing, T m j x < t c i s , nor a crossover at r m ; x ~ t c i s . However a 
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very vast experimental data set is available for a somewhat 
different system comprising a methanol solution of an am- 
phiphilic diblock copolymer (polybutylacrylate-b-polyacrylic 
acid, each block 7.5 kg/mol) quenched in water. By using a 
highly efficient impinging jet mixer, FigureQJb), Johnson and 
Prud'homme |01 succeeded in covering a very broad range 
of mixing times and observed all the three above-mentioned 
regimes. Their original data - the hydrodynamic diameters 
of the micelles formed vs the mixing time - are shown in the 
inset of Figure [3] The coagulation in a dispersion containing 
diblock copolymers must obey kinetics very similar to the one 
described by ([TJ and, thus, yield scaling © for the NP size. 
This implies that a master curve must be obtained in Figure [3] 
if one shifts the data along the horizontal axis by the poly- 
mer mass fraction w p . Moreover, a typical diameter scaling 
(TmixWp) 1 ^ is expected to be observed at long mixing times. 

One important difference between the concentration depen- 
dence of the size predicted in this Letter and the measurements 
in 10] is the fact that our equation © does not take into ac- 
count the size of the surfactant layer on top of a NP. Indeed, 
such an approximation certainly holds in case of a polymeric 
surfactant. In case of diblock copolymers, however, the size 
of the hydrophilic corona surrounding the hydrophobic core 
cannot be neglected. To compute a hydrophobic core diam- 
eter from a hydrodynamic diameter of a copolymer micelle, 
we recall that the latter scales as a power 1/5 of the micelle 
mass 01511 . As the core of a micelle consists almost solely of 
the hydrophobic polymer segments, the core size scales as a 
power 1/3 of the mass, yielding R COK oc dJ 3 . The data re- 
drawn in D^ 3 vs T m [ x c p coordinates, Figure [3] indeed shows 
a master curve obeying equation ©: it is characterized by a 
typical (Tmix/Tcis) 1 / 3 scaling at long mixing times and shows 
a plateau in the fast mixing regime, exactly as the theory pre- 
dicts. 

Note, that the NP size in Figure[3]is completely determined 
by the kinetics and is not related to the equilibrium diblock 
copolymer micelle size. Indeed the latter would depend solely 
on the molar mass, composition and solvent quality, whereas 
the NP size is a strong function of concentration. Although 
the NP system is not in a thermodynamic equilibrium, it is 
long-lived. As an X-ray study on a somewhat different di- 
block copolymer system shows iflill . micellization of copoly- 
mers includes two stages. The first rapid stage is totally con- 
trolled by kinetics and leads to NP formation described in the 
present work. The second, several orders of magnitude slower 
process, drives the NP system to the thermodynamic equilib- 
rium. 

Based on the experimental evidence discussed above, one 
can conclude that the nanoprecipitation model based on a dif- 
fusion limited coalescence mechanism adequately describes 
the NP formation process. Two relevant time scales, the mix- 



ing and the encounter and coalescence times, are identified 
in (0 and their ratio is shown to be of a critical importance 
for the NP final diameter. The latter is predicted to scale in a 
universal manner and be sensitive predominantly to the mix- 
ing time and the polymer concentration if the surfactant con- 
centration is sufficiently high. The molar mass of the carrier 
polymer is shown to have little influence. Experimental data 
available corroborate the predictions of our model and provide 
a solid framework for tailoring NPs with a priori determined 
size thus avoiding a laborious trial and error approach. 

We thank Prof. R. Prud'homme for useful discussions. 
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